Hexagonal eutectic solidification patterns operating near a marginal stability point 
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We study the long-time dynamics of hexagonal directional-solidification patterns in bulk samples 
of a transparent eutectic alloy using an optical method which permits real-time observation of the 
growth front. A slow dilatation of the patterns due to a slight curvature of the isotherms drives the 
system into a permanent regime, close to the threshold for the rod splitting instability. Thus an 
apparently minor instrumental imperfection suffices to maintain the system near a marginal stability 
point. This answers the long-standing question of spacing selection in bulk eutectic growth. 
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Directional solidification of nonfaceted binary eutec- 
tics {i.e. binary alloys, which are two-phased in the 
solid state, both solid phases displaying nonfaceted 
melt growth) gives rise to diffusion-controlled out-of- 
equilibrium patterns consisting of more or less periodic 
arrangements of the two solid phases. These patterns 
have spacing values A lying usually in the 1-10/im range 
for a solidification rate V of l/ims~^, and thus contain 
a large number of repeat units in bulk samples. What is 
their behavior at long solidification times -in other words, 
what are their permanent regimes, if any, at given con- 
trol parameters (V, and the alloy composition C)- is a 
long-standing open question [ll,!!!,!!,!^- Discussions have 
focussed on two apparently unrelated sets of results. On 
the one hand, theory 0, 0] supported by experiments 
in thin samples of transparent eutectics ^, lZ] indicates 
that the permanent regime of an ideal system must be a 
perfectly periodic pattern, the spacing of which depends 
on initial conditions. On the other hand, numerous ex- 
periments in bulk metallic eutectics displayed permanent 
regimes characterized by a broad spacing distribution in- 
dependent of initial conditions. The mean value A of 
this distribution varies with V approximately as V~^-^ 
This gap between theory and experiments 
is called the spacing-selection problem in bulk eutectic 
growth. 

To advance this, and other, fundamental problems in 
eutectic pattern formation, we built an optical device, 
which permits real-time observation of transparent eu- 
tectic growth fronts in bulk samples. A detailed pre- 
sentation of this device can be found elsewhere [llj. In 
brief, we use glass crucibles made of two glass plates 
separated by 0.4mm-thick plastic spacers. The cross- 
section of the inner volume, and thus the area occupied 
by growth front, is a 0.4 x 6mm^ rectangle. The crucibles 
are filled with a transparent eutectic alloy, and placed be- 
tween a cold oven and a hot oven. The thermal gradient 
(G = 8 ± 1 Kmm~^) is established by thermal diffusion 
along the sample. The samples are pulled at a rate V to- 
ward the cold side of the thermal gradient. The growth 
front lies in the 5-mm wide gap between the ovens and 



is observed through the liquid and a sample wall with 
a long-distance microscope. The microscope and the di- 
rection of lighting are tilted in the plane yz^ where z 
is the growth axis and y is the normal to the sample 
walls, and make angles of about 40° with z. The result- 
ing one-directional compression of the image is corrected 
numerically. To visualize the growth front itself (not the 
underlying solid), a dark- field image is formed with the 
light emerging from the interfaces between the liquid and 
one of the eutectic solid phases. Two additional refine- 
ments in the preparation of the experiments are impor- 
tant. Single eutectic grains are grown to eliminate the 
perturbations that are generated by the anisotropy of the 
surface tensions, in particular, near eutectic grain bound- 
aries jUj. The thermal bias G^y? (^^^ wanted) 
^-component of the thermal gradient, is made as small as 
possible in the region of the growth front using a device 
included in the ovens for this purpose. 

This method of observation was successfully used dur- 
ing a study of lamellar (banded) eutectic growth pat- 
terns [13]. In this Letter, we employ it to study the 
spacing-selection problem in the rod-like (hexagonal) 
growth patterns of the alloy succinonitrile-(D) camphor 
(SCN-DC) at eutectic concentration (14 mol% of DC). 
The basic properties of this alloy were studied experi- 
mentally [ij]. The scaling constant for eutectic growth 
(A^V = 10.2 ± 1.5//m^s-\ where A ra is the minimum- 
undercooling spacing [1]) was determined, and it was 
shown that anisotropy effects are particularly weak in 
this alloy. The eutectic solid phases are almost pure 
cubic SON and hexagonal DC, respectively. The en- 
velope of the growth front nearly follows the isotherm 
corresponding to the eutectic temperature of the sys- 
tem {Te = 38.3°C). The growth patterns are arrays 
of disks of DC-liquid interfaces ("DC caps") embedded 
in a continuous SCN-liquid interface. We selected the 
light emerging from the DC caps to form the images of 
the growth patterns, which thus appeared as arrays of 
bright spots on a dark background (Fig. [T]). These spots 
correspond to the caustics of the light rays transmitted 
through the DC rods and give the position of the centers 
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of the DC cap, but not their contours. The effective res- 
olution of this method of observation is of about 3/im, 
which imposes us to use sohdification rates in the 0.01- 
0.1/ims~^range in the case of SCN-DC. 




FIG. 1: a) Dark- field image of a growth front of a eutectic 
SCN-DC alloy directionally solidified at V = 0.035/ims"\ 
The growth direction z is pointing towards the reader. The 
normal y to the sample walls is pointing upwards. The ver- 
tical side of the field of view coincides with the 400/im-long 
width of the growth front, b) Voronoi diagram of the growth 
pattern. Cells have been greyed according their nearest- 
neighbor numbers (medium grey: 6- fold coordinated cells). 




FIG. 2: Left: Superposition of two binarized micrographs 
of a growth pattern at times ti (black dots) and ti + 30min 
(white dots). F = 0.035/ims~^. Vertical dimension: 400/im. 
Right: Sketch of a longitudinal section of the system assuming 
a curved growth front. 

We registered the trajectories of all the DC caps ap- 
pearing at the growth front during long solidification 
runs. This allowed us to reveal that, in all the exper- 
iments, the DC caps were drifting slowly along the nor- 



mal to the sample walls (Fig. [2]). No drift parallel to the 
walls was observed. The drift velocity Vd was a definite, 
roughly linear function of y that went to zero at a position 
which was slightly sample dependent, but was always 
located inside the sample. Thus the spatiotemporal di- 
agram of the patterns was fan-shaped, and the patterns 
were undergoing a permanent dilatation along y. A sim- 
ple explanation for this observation is that, in our setup, 
the isotherm at Te has a negative curvature (it is bulging 
toward the liquid) in the plane yz and that the trajecto- 
ries of the DC caps remain locally perpendicular to this 
isotherm (Fig. Ob). While the latter statement, known 
as the normal-growth condition, is certainly correct to a 
good approximation, the former (a negative curvature of 
the isotherm) needs further discussion. 

The curvature of the isotherms in our system is slight, 
as will be seen shortly, and could not be measured by di- 
rect observation. We determined it from the drift of the 
growth pattern as follows. According to our conjecture, 
Vd/V = dC/dy^ where z = ({y) is the equation of the Te 
isotherm in the plane yz. A linear fit on the measured 
values of Vd{y) gave us a linear function for d^/dy^ and, 
by integration, a parabolic function for Civ)- Tit turned 
out that the radius of curvature R at the vertex of the 
parabola {i.e. at = yo) was always much larger than the 
sample thickness. The curvature was thus independent 
of y within the experimental uncertainty. Measurements 
of R performed with this method in nine samples for dif- 
ferent values of V and A ranged from about 2.2 to 5.4mm. 
Incidentally, an elementary calculation gives the residual 
thermal bias knowing R^ yo and G. The maximum value 
of Gy was about 0.4Kmm~^ in our experiments. 

We obtained another estimate of R from the rate of in- 
crease of A over time. We define A and A as the mean val- 
ues of the nearest-neighbor distances for a given DC cap 
and the whole growth front, respectively, at a given time 
t. An elementary geometrical calculation shows that, for 
an isotherm of uniform one-directional curvature R~^, 
A ^ Aoe^/'^, where Aq is the value of A at t = and 
r = 2R/V. This equation is valid when the number 
of repeat units is locally conserved during the process, 
i.e. when no DC cap is created or eliminated. Whether 
or not this conservation condition was fulfilled during a 
given portion of a solidification run could be checked by 
direct observation. DC cap elimination events were never 
observed (except, of course, at the glass walls) during this 
study, as could be expected since the pattern was under 
stretching, and can be left aside. Figure [3] shows the time 
evolution of A during a solidification run without creation 
events. An exponential fit (which was practically linear 
in the time interval of interest) yielded r = 82 ± 12h, 
and hence R = 3.2 ± 0.5mm, in good agreement with the 
value of R derived from the measurement of Vd{y) during 
the same experiment. 

What is the origin of the curvature of the isotherms 
in our directional-solidification setup? Several factors 
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could be at play, among which the differences in thermal 
conductivity of the various media composing the sam- 
ples, the rejection of latent heat at the growth front, and 
convection-induced transverse gradients. However, only 
the first of these factors (non-uniform thermal conduc- 
tivity) is likely to be important in our case because of 
the very low solidification rates used, and the order of 
magnitude of ratio of R to the sample thickness (0.1) ob- 
served. In any case, it should be noted that the effect 
of a slight curvature of the isotherm on growth pattern- 
ing presumably reduces to the normal-growth condition, 
which is independent of the origin of the curvature. 
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FIG. 3: Time variation of the mean spacing in the absence 
of rod splitting events. V — 0.021/ims~^. Open circles: mea- 
sured A values. Continuous line: exponential fit with a char- 
acteristic time of 82/i. Insets: Voronoi diagrams of a fixed 
region of the growth pattern at the beginning (upper left cor- 
ner) and the end (lower right corner) of the process. 

Bifurcations to oscillatory regimes apart, the upper 
stability bound of hexagonal eutectic patterns at fixed 
V corresponds to an instability leading to the creation of 
new repeat units, and hence to a decrease in A [l5[. This 
instability was observed in our experiments at sufficiently 
long solidification times. It consists of the splitting of 
a DC cap into two parts (Fig. [4j). A detailed report 
on this, and other instability processes (rod termination, 
oscillatory bifurcations), will be presented in the future. 
For our present purpose, the important point is that, 
in long-duration experiments at constant V ^ the average 
rod splitting frequency eventually adjusted itself so as 
to counterbalance the increase of A driven by the cur- 
vature of the isotherms. Any solidification rate program 
ending with a maintain at constant V led to a permanent 
regime characterized by a plateau in the \{t) curve. Illus- 
trative examples are shown in Figure [5l Direct counting 
of the rod splitting events confirmed that this plateau 
corresponded to a global balance between rod splitting 
and the curvature-driven increase of A. In this regime, 
A fluctuated with characteristic times on the order of an 
hour about a constant value \p. Data from seven exper- 
iments performed with different solidification rate pro- 
grams yielded A^ ex y-^/^ and A^/A^ = 1.03 ± 0.04. A 
quantitative comparison between this quantity and the 



threshold spacing for the rod splitting instability is be- 
yond the scope of this article. However, it is worth noting 
that the local spacing values at the onset of rod splitting 
events were slightly larger than Ap. 




FIG. 4: Direct observation of rod-splitting events. Snapshots 
of a 65 X 74/im^ region of a growth pattern taken at time 
intervals of 7 min. V = 0.035/ims-\ Two DC caps (located 
in the second column from left) are simultaneously splitting. 

To sum up, long-duration directional-solidification 
runs at constant V comprised two stages: first, a tran- 
sient, during which the mean spacing A increased or de- 
creased depending on its initial value; second, a perma- 
nent regime, in which A fluctuated about a constant value 
\p slightly smaller than the threshold spacing for rod 
splitting. What was the evolution of the degree of or- 
der of the growth patterns during these stages? During 
transients without rod splitting (Fig. [5la), we observed a 
continuous ordering process, consisting in the formation 
and growth of perfectly ordered domains separated by 
thin boundaries (Figs. [T]and[3]). Decelerating the solidifi- 
cation during such a transient made it last longer, which 
furthered the ordering process [16]. The local mecha- 
nisms of the ordering process remain to be studied. We 
think that the curvature-driven stretching of the pattern 
played an important role in this process, in particular, 
because the hexagonal array had a nearest-neighbor di- 
rection either parallel or normal to the sample walls in- 
side most domains. The order that was established dur- 
ing transients without rod splitting was progressively de- 
stroyed when rod splitting set in. So long as there re- 
mained distinct domain boundaries, rod splitting events 
occurred preferentially in these boundaries, and had some 
short-range correlation. The more or less periodic fluc- 
tuations of A observed in the well-developed permanent 
regime are perhaps due to collective phenomena (such as 
the "cascades" observed in curved cellular growth fronts 
[13]), but we could not establish this clearly. For com- 
parison with metallurgical studies, the time evolution of 
the A distribution during a representative solidification 
run is shown in Figure [6l To note, the A distribution be- 
comed peaked about Am, but still extended over a wide 
range in the permanent regime. 

In conclusion, we have observed, in a transparent 
eutectic directional-solidification system, a permanent 
regime, which presents statistical features (sharp selec- 
tion of the mean spacing, broad dispersion of the lo- 
cal spacing values) similar to those that were previously 
noted, but not understood, in metallic eutectics. We have 
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FIG. 5: Time variation of the mean spacing A during two dif- 
ferent solidification rate programs. The time variations of V 
and A are shown in the upper and lower parts of the diagrams, 
respectively. The initial value of the mean spacing belonged 
to the stability interval in (a), and was above this interval in 
(b). In the latter case, rod splitting occurred from the start, 
and became less and less frequent until a permanent regime 
was settled. 
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FIG. 6: Normalized histogram of the local spacing at the 
indicated times during the solidification run of Fig. Xm- 
minimum-undercooling spacing 



brought to light the dynamical origin of this permanent 
regime. A slight instrumental imperfection -namely, a 
negative curvature of the isotherms- drives the system 
towards the marginal stability point for the repeat-unit 
creation instability of the growth pattern. After a tran- 
sient, the two processes (external forcing and dynamical 
instability) come into balance, and the spacing remains 
close to the instability threshold value. In general terms, 
this conclusion is not a surprise. Langer [2] pointed out 
long ago that the broad multistabilty of ideal directional- 
solidification systems could be easily lifted by some exter- 
nal noise or forcing, and that, then, the system was likely 
to operate near a marginal stability point. Our study per- 
mits to specify that a curvature of the isotherms was the 
external forcing at play in our system. A similar forcing 
is certainly at play in all nonfaceted eutectic directional- 
solidification systems, but it may be combined with other 
factors capable of producing a drift of the pattern such as 
capillary anisotropy and a transverse thermal bias. These 
factors were both small in our experiments, as mentioned, 
which allowed us to observe the effect of a curvature of 
the isotherms in their plainest form. 
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